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HYDROFOIL  SHIP 

CONSTANT  LIFT  CONTROL  SYSTEM 

PART  I:   HYDRODYNAMICS 

MICHAEL  ROY  TERRY 

Submitted  to  the  Departments  of 
NAVAL  ARCHITECTURE  and  MARINE  ENGINEERING 
and  MECHANICAL  ENGINEERING  on 
23  May  1969,  in  partial 
fulfillment  of  the  requirements 
for  the  Degrees  of  NAVAL  ENGINEER, 
and  MASTER  of  SCIENCE  in  MECHANICAL  ENGINEERING. 


This  paper  analyzes  the  unsteady  hydrodynamics  of 
1)  FULL  INCIDENCE;  2)  TRAILING  EDGE  FLAP;  and  3)  TAB  FOIL 
hydrodynamic  control  devices.  The  device  is  to  be  used  in  a 
constant  lift,  wave  alleviation  control  system  on  an  OPEN- 
OCEAN  HYDROFOIL  SHIP.   Each  device  is  compared  in  POWER,  DRAG, 
and  WEIGHT.  The  conclusions  reached  are: 

1.  Wave  Alleviation  is  an  important  aspect  of 
hydrofoil  ship  control  system  design. 

2.  Unsteady  Analysis  is  vital  to  understanding 
the  response  of  a  foil  in  waves. 

3.  Negative  Control  System  Power  is  required  for 
wave  alleviation  by  a  tab  fol*  with  a  lagging 
tab. 

4.  Trailing  Edge  Flap  requires  less  overall 
system  power  thah  a  Full  Incidence  Foil. 

5.  Hydrofoil  Ship  Operation  Policy  will  be  greatly 
affected  by  unsteady  hydrodynamics. 

6.  Unsteady  Foil  Response  is  not  a  unique  function 
of  encounter  frequency. 
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SYMBOLS 

a  ■  NACA  camber  loading  parameter 

a  ■  a  hinge  point  in  %   of  b 

A  ■  Foil  area 

b  *  Semichord 

c  ■  Chord 

c  =3  hinge  point  in  %   of  b 

C  *  Designates  circulatory 

C-l  -  Lift  coefficient   L/( (l/2)pV2 A) 

C,0  *  Lift  coefficient  at  ideal  angle  of  attack 

C(k)  =  Theodorsen  2-D  unsteady  lift  function 

d  «  Source  position 

D  ■  Drag 

e  =  Designates  encounter 

g  ■  Gravitation  constant  ft2/sec 

h  e  Keel  to  foilborne  waterllne  height 

1  m   /=X 

I_  «•  Main  foil  inertia/ ft  span 

Ib  ■  Flap  inertia/ft  span 

Jn(k)  ■  Bessel  Function  of  order  n 

k  ■  Reduced  frequency 

L  *  Li ft /ft  span 

L  »  Length 

m  *  Source  strength/ft  span 

M  ■  Moment /ft  span 

M  ■  Mass  of  the  ship 


NC  ■  Denotes  non-circulatory 

Q  ■  Flow  rate 

s  ■  Span 

t  ■  Time 

t  =  Thickness 

T  =  Period 

U  -  Ship  Speed 

V  =  Wave  Phase  Velocity 

w  ■  Denotes  wave 

x  ■  Distance  in  horizontal  direction 

x0  =  Longitudinal  distance  from  reference  point  to 
center  of  gravity 

y  ■  Distance  in  vertical  direction 

z  ■  Distance  traveled  in  vertical  direction 

z  =  Vertical  velocity 

M 

z  *  Vertical  acceleration 

Z  ■  Vertical  force 

a  =  Pitch  angle  of  incidence  about  ba 

a0  -  Pitch  angle  of  incidence  magnitude 

a  «  Incidence  angular  velocity 

a  ■  Incidence  angular  acceleration 

0  ■  Pitch  angle  of  flap  about  be 

30  «  Pitch  angle  of  flap  magnitude 

• 

3  ■  Flap  angular  velocity 

3  ■  Flap  angular  acceleration 
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A  ■  Ship  displacement 

A  ■  Increment  of  change 

e  «=  Small  increment 

nD  =  Propulsive  coefficient 

n  ■  Wave  height 

Ho  e  Wave  height  magnitude 

r\  ■  Wave  orbital   velocity 

n  ■  Wave  orbital  acceleration 

M 

9  *  Ship  pitch  acceleration 

X  =  Wave  length 

tt  =  3. 1^16... 

p  -  Fluid  density  in  slugs/ft8 

<fr  ■  phase  angle 

A<f>  ■  Phase  angle  between  a,  3 

ty  =  Angle  of  attack  of  ship  to  wave3 

^  a  Stream  function 

a)  «  Circular  frequency  radians/sec 
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MAJOR  COMPUTER  PROGRAM  SYMBOLS 


Plotting  symbol  for  |a| 


ALPHA 

s 
a 

a 

a 

ALP HA DOT 

a 

ALPHADDOT 

a 

B 

a 

b 

B 

a 

Plotting  symbol  for  |$| 

BA 

SI 

s= 

3 

s 

30/a0 

BETA 

B 

BETADOT 

BETADDOT 

3 

C 

a 

c  hinge  line  in  %   of  b 

CA 

a 
a 

k 

CK 

C(k) 

DPHE 

a 

A<J> 

PI  A 

s 

Xa 

FIB 

a 

Jb 

HW 

s 
a 

n 

LIPTW 

Lw 

LI  TO  A 

L 
a 

LIPTB 

H 

LIPTOT 

MOMTW 

KOMTBW 

MOMTA 


■  Total  lift  per  unit  span 

■  \\j   about  ba  due  to  waves 

■  Mw  about  be  due  to  waves 

■  Ma  about  ba  due  to  a 
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riOMTbA  ■  W a  about  be  due  to  a 

M0MT13  ■  Mg  about  ba  due  to  3 

MOMTBB  -  M&  about  be  due  to  3 

MOMTOT  ■  Total  hydrodynamic  moment  about  ba  per  unit 
span 

riOMTBT  ■  Total  hydrodynamic  moment  about  be  per  unit 
span 

MOMAPA  ■  Required  applied  moment  about  ba  per  unit  span 

MOMAPB  »  Required  applied  moment  about  be  per  unit  span 

P  »  Plotting  symbol  for  PAVEA 

PAVEA  «  Average  power  about  ba  per  unit  span 

PAVEB  «  Average  power  about  be  per  unit  span 

PHEW  =  *w 

Q  »  Plotting  symbol  for  PAVEB 

RO  a  p 

31  =  y 

UPW  a  Upwash  velocity   r\ 

V/E  a  0) 

WL  a  X 

Y  a  Plotting  symbol  for  ^3 

Z  a  Plotting  symbol  for  ^a 

A  =  a  hinge  line  in  %   of  b 
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I.   INTRODUCTION 

Renewed  interest  in  the  use  of  hydrofoil  lifting  sur- 
faces to  raise  a  ship  above  the  surface  of  the  water 
occurred  in  the  1950* s.   Since  that  time,  development  of 
complete  hydrofoil  ship  systems  has  progressed  at  a  rapid 
rate.   During  the  recent  development  stages,  the  encounter 
between  ships  and  open-ocean  vraves  has  identified  the  pro- 
blem of  ship  motion  control  as  being  critical.   The  foil 
encounters  lift  changes  caused  by  wave  particle  orbital 
motion.   The  wave  particle  orbital  velocity  is  superposed 
on  the  uniform  in-flow  velocity  (ship's  speed  through  the 
water).   The  superposition  results  in  an  angle  of  attack 
change  relative  to  the  foil  and,  therefore,  changes  its 
lift.   The  changes  in  lift  are  disturbing  forces  to  the 
ship  and  occur  at  the  encounter  frequency  of  the  ship  meeting 
the  waves.   Figure  (1)  depicts  this  situation. 
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crest  to  trough 


keel  to  foil 
V 


figure    1  7 

This  paper  will  address  only  open-ocean  waves  and  ocean- 
going hydrofoil  ships. 

During  the  early  stages  of  development,  attention  was 
focused  on  the  choice  of  foil  configuration  for  best  results. 
From  that  lively  discussion  and  experimentation,  two  major 
design  configurations  developed.   The  first  is  commonly 
referred  to  as  surface  piercing,  and  the  second  as  f ul ly 
submerged.   Figure  (2)  depicts  these  configurations. 


surface    piercing 


fully  submerged 
figure     2 


15 


This  paper  will  address  only  the  fully  submerged  configura- 
tion. Hydrofoil  ships  with  fully  submerged  foils  are  inher- 
ently unstable  in  heave,  pitch,  and  roll  when  the  hull  is 
out  of  the  water,  because  the  foil  has  very  little  sensiti- 
vity to  depth  below  the  water  surface.   Therefore,  a  control 
system  must  be  provided. 

While  foilborne,  two  modes  of  operation  are  possible. 
The  first  is  the  platform  mode  where  the  ship's  motion  is 
maintained  at  constant  height  above  the  mean  water  level, 
at  constant  attitude  about  the  height,  and  with  no  accelera- 
tions.  The  platform  mode  is  only  possible  in  wave  heights 
that  do  not  exceed  the  keel-to-foilborne  waterline  distance. 
The  second  is  the  contour  mode  where  the  ship's  motion  is 
maintained  at  constant  height  above  the  local  water  surface 
at  two  points  along  the  ship's  length.   The  contour  mode 
is  only  necessary  in  wave  heights  that  exceed  the  keel-to- 
foilborne  waterline  distance.   Figure  (3)  depicts  these 
modes . 


contouring 


figure     3 
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platforming 


Thia  paper  will  address  only  the  platform  mode  and  the  wave 
systems  that  allow  the  use  of  that  mode. 

The  design  of  motion  control  for  a  hydrofoil  ship  must 
meet  the  three  basic  objectives  outlined  In  Figure  (4). 


STABILITY 


MANEUVERABILITY 


WAVE  ALLEVIATION 


(Maintaining  Commanded 
Equilibrium  Attitude) 


(Changing  to  a  New 
Commanded  Equilibrium) 


(Reducing  the  Effect 
on  Lift  due  to  Waves) 


figure      4 


Most  attention  in  the  past  has  been  concentrated  on  stability 
and  maneuverability  control  functions.   The  effect  of  waves 
on  the  ship  has  been  treated  as  a  disturbance  from  equili- 
brium, and  the  stability  elements  of  the  control  system  have 
been  sized  to  be  able  to  correct  the  ship's  attitude  and 
reduce  the  accelerations  to  a  tolerable  level.   In  contrast, 
this  paper  will  address  only  the  wave  alleviation  objective. 

The  control  system  is  usually  composed  of  elements  as 
shown  in  Figure  (5). 
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SENSORS 


(ELECTRONIC) 


AUTOPILOT/COMPUTER 


(ELECTRONIC) 


SERVO 


(FLUID  POWER) 


HYDRODYNAMIC  DEVICE 


(LIFTING  SURFACE) 


figure       5 


This  paper  will  address  only  the  Hydrodynamic  Device, 


13 


II.   PROJECT 


The  objective  of  the  above  focus  is  now  described.   An 
open-ocean  hydrofoil  ship  will  operate,  at  times,  when  wave 
systems  are  present  for  the  complete  voyage.   Operation  in 
waves  is,  therefore,  an  important  aspect  of  the  design.   Would 
the  control  system  perform  better  than  present  day  systems 
if  the  approach  to  its  design  considered  only  wave  allevia- 
tion as  the  objective?  The  hypothesis  i3  that  the  power 
level  of  the  control  system  operating  in  waves  could  be 
significantly  reduced  if  this  approach  was  taken.   Also,  the 
wave  alleviation  control  system  could  be  integrated  with  a 
present  day  standard  stability  and  maneuverability  control 
system  to  produce  a  control  system  with  overall  performance 
superior  to  today's  design.   The  stabilizing  and  maneuvering 
elements  would  effectively  be  operating  as  if  they  were  in 
calm  water,  when  in  fact  the  hydrofoil  ship  is  flying  in 
open-ocean  waves,  Figure  6. 


A 


ift 


inflow 


t> 


J 


figure     6 
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In  the  process  of  designing  the  optimum  control  system  for 
wave  alleviation,  different  solutions  must  be  tested.   The 
criteria  best  suited  for  this  test  is  minimum  overall  system 
power.  The  three  elements  that  contribute  directly  to  over- 
all system  power  are  shown  in  Figure  (7). 


POWER 


Power  to  Operate 
the  Control  Surface 


DRAG 


Power  to  Overcome  the  Drag 
of  the  Hydrodynamic  Device 


WEIGHT 


Power  to  Lift,  Through  the 
Lift/Drag  Ratio,  the  Weight 
of  the  Control  System 


figure     7 


For  different  designs,  any  one  of  the  above  elements  of 
overall  control  system  power  may  predominate.   The  task  of 
designing  the  control  system  can  be  broken  down  into  the 
design  of  the  four  basic  elements — sensor,  autopilot,  servo, 
and  hydrodynamic  device  as  shown  in  Figure  (5)»   "Hydrofoil 
Ship  Constant-Lift  Control  System:   Part  I  Hydrodynamics" 
will  address  only  the  fourth  element— Hydrodynamic  Force 
Producing  Device. 
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In  summary,  Figure  (8)  shows  the  steps  in  focusing  this 
thesis  within  the  overall  project. 
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figure    8 
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III.   IIYDRODYNAMIC  DEVICE 
A.   DEVICES 

The  first  step  in  the  plan  to  find  the  minimum 
power  device  is  to  develop  as  many  candidates  as  possi- 
ble for  closer  scrutiny.   In  general  terms,  the  hydro- 
dynamic  device  must  effect  a  force.   Although  there 
are  a  number  of  ways  of  effecting  a  force,  this  paper 
will  consider  only  those  that  develop  lift.   This  is  a 
reasonable  specialization  because  the  hydrofoil  ship 
utilizes  lifting  surfaces  for  lifting  itself  out  of  the 
water  and  for  stability  and  maneuverability. 

Table  (1)  is  an  attempt  at  classifying  the  candi- 
dates .   This  paper  will  pursue  only  the  following 
candidates: 

1.  Full  Incidence 

2.  Trailing  Edge  Flap 

3.  Tab  Foil 

Appendix  (A)  describes  the  remaining  candidates  and 
discusses  some  of  their  characteristics  in  qualitative 
terms. 
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B,   SHIP 

before  proceeding  v/ith  a  quantitative  analysis  of  the 
three  candidates,  a  reference  hydrofoil  ship  3hould  be  de- 
fined for  use  as  a  common  denominator  for  each  design. 
Table  (2)  lists  the  parameters  of  the  reference  hydrofoil 
ship  chosen. 

TABLE  (2) 

Displacement 

Length  Overall 

Design  Speed 

Design  Lift/Drag 

Design  Propulsive  Coefficient  n  «  .6 

P 

Each  of  2  Main  Lifting  Foils  NACA  16-408 
Thickness/Chord         t/c  ■  .08 
Mean  Line  a  ■  1.0 

Design  Lift  Coefficient  Cjo  B    .4 
Chord  c  =  8.58  ft 

Semichord  b  =  4.29  ft 

Span  s  «  21.25  ft 

Weight  wt.  ■  7  tons 

Foil  Inertia/Ft  Span    la   =  140.0  Slug  Ft2/ft 

1/4  Chord  Trailing  Edge  Flap 

Chord  =  2.145  ft 

Wt.  »  1.75  tons 

Flap  Inertia/Ft  Span    lb   =8.80  Slug  Ft2 /ft 

Keel  to  Foilborne  Waterline  Height  h  e  14.0  ft 


IV.   POWER 

The  first  and  most  Important  component  of  overall  con- 
trol system  power  Is  the  power  to  drive  the  hydrodynamlc 
device  through  motions  that  will  effect  platformin^  in 
waves.   A  reasonable  model  of  motion  in  the  vertical  plane 
for  the  hydrofoil  ship  is 

M 

Mz  +  Mx08  ■  Z (waves)  +  Z(Ship    +  Z( Control     (1) 

Motion)     Motion) 

The  effects  of  waves,  ship  motion,  and  control  motion  are 
modeled  as  linearly  superimposed.   For  the  platforming  mode, 
the  equation  Is  reduced  to 

Zero  ■  Z( waves)  +  Z( control  (2) 

motion) 

The  objective  is  to  cancel  the  lift  caused  by  waves  by 

generating  an  equal  and  opposite  lift  caused  by  the  control 

surface.   From  the  required  motion  of  the  control  surface 

and  the  resultant  forces  on  It,  the  power  required  can  be 

calculated 


Power  »  (motion  velocity) (summation  of  forces  in   (3) 

the  direction  of  the  velocity) 


The  next  task  is  to  look  at  the  effect  of  waves  on  a  foil. 
A.   QUASI-STEADY  LIFT  DUE  TO  WAVES 

From  a  qualitative,  quasl-3teady  point  of  view,  the 
lift  due  to  a  wave  passing  a  fully-submerged  hydrofoil  Is 
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caused  by  the  orbital  velocity  of  the  water  particles  super- 
posed with  the  velocity  of  the  foil  through  the  water.   For 
a  steady  foil,  the  orbital  velocity  appears  as  an  ancle  of 
attack  change  for  the  incoming  flow,  Figure  (9). 


A 


u 


(3;\ 


ft 


V    •^-^Oe/ 


figure    9 


For  the  hydrofoil  ship,  trie  limiting  wave  condition  for 
continued  platforming  is-- total  amplitude  of  the  wave  equal 
to  the  keel-to-foilborne  watorline  height.   For  the  parti- 
cular hydrofoil  ship  chosen,  this  means  a  wave  half-height 
of  7.C  feet.   For  an  open-ocean  wave,  the  common  wave  length 
for  this  amplitude  is  280.0  feet.   At  the  ocean  surface,  at 
a  given  point,  a  reasonable  model  of  the  wave  vertical  posi- 
tion Is 


n   -   Re{n0exp(-i2irV/X-t)> 


For  a  deep  water  wave 


V.-  /gT72T 
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or 

n  -   Re{n0exp(-l/21rc7r  •    t)>  (6) 

The  wave   vertical   particle   velocity   is   then 

fi  -  Re{n0(-i/^G7T)exp(-i/7l?g7T  •   t)>  (7) 

The  amplitude  of  the  orbital  velocity  is  then 

| n |  -  no/27g7X  (8) 

For  the  chosen  wave 

|n |   =5.85  ft/sec  (9) 

With  the  hydrofoil  ship's  velocity  at  67.5  ft/sec,  then  the 
maximum  angle  of  attack  change  due  to  vertical  orbital  velo- 
city is 

Tan  a  =  5.95/67.5  (10) 

or 

a  *  5.060  (11) 

The  maximum  inflow  velocity  change  due  to  horizontal  orbital 
velocity  is 

Ui  =  67.5  ±  5.95  ft/sec  (12) 

To  determine  the  effect  on  lift,  we  must  look  at  the  lift 
coefficient  C^,      For  the  5.06°  angle  of  attack  change,  flat 
plate  model  of  the  foil  gives 
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3(^/30  *  2ir  (13) 

or 

AC^  ■  2irAa  (1M) 

ACX  «  .5'46  (15) 

so 

c!i  e  .^  +  .5^6  »  .9^6  (16) 

For  the  hydrofoil  ship  chosen,  the  ratio  of  the  lifts  is  then 

I^/L,  =  2.36  (17) 

In  contrast,  for  the  5.95  ft /sec  inflow  velocity  change ,  the 
ratio  of  the  lifts  is 

L4/L0  =»  (Uj/U0)2  (18) 

Lj/L0  *  1.19  (19) 

Prom  the  above,  the  effect  on  lift  due  to  the  vertical  compo- 
nent of  orbital  velocity  is  far  greater  than  that  due  to  the 
horizontal  component  of  orbital  velocity.   This  paper  will 
only  consider  the  effect  of  the  vertical  component  of  orbital 
velocity. 

The  use  of  quasi-steady  analysis  is  not  adequate  for 
the  case  of  the  hydrofoil  ship  chosen.   The  encounter 
frequency,  semichord,  and  ship's  velocity  combine  to  give 
a  reduced  frequency 

k  =web/U  (20) 
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in  the  range  where  unsteady  effects  are  most  important, 
as  will  be  shown  in  Section  (B).  Figure  (10)  describes 
the  concept  of  encounter  frequency  u>e. 


here 


figure    10 


we  =    2ir/Te   ■   27iUe/Le 


}/Cos(fJ 


(21) 


where 


and 


Le    a    X/Cos\j/ 


(22) 


U0   =    (V/Cosi|»)    +    U 


(23) 


so 


u)e   -    (2ir/X)(V  +    UCosi};) 


(21) 


but 


V2    =    gX/2ir 


(5) 


therefore 


X   -   2ng/u)w2 


(26) 
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or 


w«  **  w  ,  +  tu  2  •  UCos(j;/g 

6      W      w 


(27) 


For  the  wave  and  ship  chosen 


we  «  .85  +  .72(67.5)(D/32.2  =  2.36  Rad/Sec 


(28) 


Semi chord     b   =   M . 29   ft 
Speed  U  =   67.5   ft/3ec 


so 


k  *   2.36(4.29)/67.5  =    .15 


(29) 


or  k  is  in  the  most  critical  range  for  unsteady  effects  as 
shown  in  Section  B. 


B.   UNSTEADY  LIFT  DUE  TO  V/AVES 


U 


MoMrBW 


figure  11 


31 


The  following  ±3  an  adaptation  of  Chapter  5  of 
Reference  (1)  for  this  particular  problem.   Please  refer 
to  Table  (3)  in  this  paper  for  the  following  discussion. 
A  model  of  a  wave  passing  a  steady  foil,  shown  in  Figure  (11), 
i3  expressed  by  Equation  (30).   The  disturbing  upwash  velo- 
city UPW  is  equal  to  a  complex  phasor  with  magnitude — |  n  | , 
phase  angle  (relative  to  wave  height)  —  <J>W,  and  being  driven 
around  the  phasor  diagram  with  angular  velocity— u)Q. 

Assuming  negligible  effects  from  free  surface,  negli- 
gible energy  radiation  on  the  free  surface;  idealized  2-D 
flat  plate  theory,  with  no  separation  or  cavitation,  the 
unsteady  lift  is  expressed  as  the  quasi-steady  lift  modified 
by  the  Theodorsen  Function  and  additional  Bessel  Functions. 
These  functions  are  expressed  in  Equations  (31  to  3*0.   The 
Theodorsen  Function  is  shown  in  Figure  (12). 

5 


1.0 


Im 


figure   12 


From  Figure  (12),  it  can  be  3een  that  the  Theodorsen  Function 
reduces  in  magnitude  (from  1  to  1/2)  and  has  its  maximum 
lagging  phase  angle  change  near  k  ■  .15; — exactly  where  the 
hydrofoil  ship  is  operating. 
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The  expression  for  the  complex  phasor  LI FTW  is  shown 
in  Equation  (35).   The  complex  phasor  moment  MOMTW  due  to 
LI FTW  acting  about  point  ba  (Pi cure  (11))  is  expressed  in 
Equation  (36).   The  complex  phasor  moment  MOMTBW  due  to  lift 
on  the  flap  acting  about  point  be  (Figure  (11))  is  expressed 
in  Equation  (37).   These  functions  are  best  visualized  with 
a  phasor  diagram.   Figure  (13)  shows  a  phasor  diagram  scaled 


to  illustrate  the  relationship  between  UPW,  LIFTW,  MOMTW, 


and  MOMTBW  for  u>e  =  (  +  ). 
NOTE : 

1.  The  diagram  rotates  clockwise. 

2.  The  unwa3h  velocity  leads  the  wave  height  by  90°. 


3.  LIFTW  lags  UPW. 

4.  MOMTW  is  in  phase  with  LIFTW. 


5.  M0MT3W  is  130°  out  of  phase  with  LIFTW. 

6.  LIFTW  and  MOMTW  would  have  been  larger  in  magnitude 
and  in  phase  with  "DW,  if  k  ■  0. 

Figure  (1*0  shows  a  phasor  diagram  again  scaled  to  illustrate 
the  same  relationships  for  w  =  (-). 
NOTE : 

1.  The  diagram  rotates  count erclockwise. 

2.  The  upwash  velocity  lags  the  wave  height  by  90°. 


3.  Again  LIFTW  and  MOMTW  lag  UPW. 

4.  MOMTBW  is  180°  out  of  phase  with  LIFTW. 
Figure  (15)  illustrates  the  meaning  of  we  ■  (+,-). 
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figure   14 


35 


m 


D 
O) 


C.   UNSTEADY  LIFT  DUE  TO  \M  OSCILLATING  FOIL 


u 


figure  16 


1.   Alpha  Motion 

The  following  is  again  an  adaptation  of  Chapter  5  of 
Reference  (1).   Please  refer  to  Table  (M)  in  this  paper  for 
the  following  discussion.   A  model  of  full  incidence  motion 
(Alpha  Motion)  (Figure  (16))  is  expressed  in  Equations  (38 
to  MO) .   It  is  expressed  as  a  complex  phasor  with  phase 
angle  relative  to  the  wave  height— $!,  and  being  driven 
around  the  phasor  diagram  with  angular  velocity — we. 

Again,  assuming  negligible  effects  from  free  surface, 
negligible  energy  radiation  on  the  free  surface;  idealized 
2-D  flat  plate  theory,  with  no  separation  or  cavitation, 
the  unsteady  complex  phasor  lift  (LiFTA)  is  expressed  in 
non-circulatory  velocity  and  acceleration  terms,  and  circu- 
latory position  and  velocity  terms  modified  by  the  Theodorsen 
Function.   LI FT A  Is  expressed  In  Equation  (Ml).   The  unsteady 
complex  phasor  moment  MOMTA  due  to  LI FT A  acting  about  point 
ba  (Figure  (16))  is  expressed  in  Equation  (42).   The  unsteady 
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complex  phasor  moment  MOMTBA  due  to  lift  on  the  flap  acting 
about  point  be  (Figure  (16))  is  expressed  in  Equation  (**3). 
Again,  these  functions  are  best  visualized  with  a  phasor 
diagram.   Figure  (17)  shows  a  phasor  diagram  scaled  to 
illustrate  the  relationship  between  ALPHA,  ALPHADOT,  ALPHA- 


DTOT,  LIFTA,  MO.MTA,  and  MOMTBA  for  we  «  (+). 
NOTE: 

1.  The  diagram  rotates  clockwise. 

2.  ALPHADOT  leads  ALPHA  by  90°. 

3.  ALPHADDOT  leads  ALPHADOT  by  90°. 


l\.      LIFTA  lags  ALPHA. 


lj.      MOMTA  and  MOMTBA  lag  LIFTA. 

6.   LIFTA  will  be  in  phase  with  ALPHA  for  k  =  0 


because  although  LIFTA  has  contributions  from 


ALPHADOT  and  ALPHADDOT,  these  are  zero  for  k  =  0, 
Figure  (18)  shows  a  phasor  diagram  to  illustrate  the  same 
relationships  with  we  =  (-) . 

NOTE : 

1.   The  diagram  rotates  counterclockwise. 


2.  ALPHADOT  leads  ALPHA  by  90°. 

3.  ALPHADDOT  leads  ALPHADOT  by  90°. 


4.   LIFTA  lags  ALPHA. 

2.   Beta  Motion 

The  following  is  again  an  adaptation  of  Chapter  5  of 
Reference  (1).   Please  refer  to  Table  (5)  in  this  paper  for 
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the  following  discussion.   A  model  of  trailing  edge  flap 
motion  (Beta  Motion)  (Figure  (16))  is  expressed  in  Equations 
(M  to  *l6) .   It  is  expressed  in  terms  of  the  complex  phasor 


ALPHA.   The  ratio  of  magnitude  between  BETA  and  ALPHA  is — 
3o/a0;  the  phase  difference  between  BETA  and  ALPHA  (+  for 
BETA  leading)  is— A<J>;  and  the  phasor  is  being  driven  around 
the  phasor  diagram  with  angular  velocity — we. 

Again,  assuming  negligible  effects  from  free  surface, 
negligible  energy  radiation  on  the  free  surface,  idealized 
2-D  flat  plate  theory;  with  no  separation  or  cavitation, 
the  unsteady  complex  phasor  lift  (LIFTB)  is  expressed  in 
non-circulatory  velocity  and  acceleration  terms,  and  circula- 
tory position  and  velocity  terms  modified  by  the  Theodorsen 
Function.   LIFTB  is  expressed  in  Equation  d7).   The  unsteady 
complex  phasor  moment  MOMTB  due  to  LIFTB  acting  about  point 
ba  (Figure  (16))  is  expressed  in  Equation  (^7).   The  unsteady 
complex  phasor  moment  MOMTBB  due  to  lift  on  the  flap  acting 
about  point  be  (Figure  (16))  is  expressed  in  Equation  (^9). 
Again,  these  functions  are  best  visualized  with  a  phasor 
diagram.  Figure  (19)  shows  a  phasor  diagram  scaled  to 
illustrate  the  relationship  between  BETA,  BETADOT,  BETADDOT, 
OTTO,  M5MTO,  and  M6MCBB  for  we  «  (  +  ). 

NOTE:  Same  as  Figure  (17). 
Figure  (20)  shows  a  phasor  diagram  to  illustrate  the  same 
relationships  with  u>e  ■  (-). 

NOTE:  Same  as  Figure  (18). 
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D.   REQUIRED  FOIL  MOTION 


Mode  1  ■  ALPHA  only  (Full  Incidence) 
Mode  2  ■  BETA  only  (Trailing  Edge  Flap) 


Mode  3  ■  ALPHA  and  BETA  (Tab  Foil) 

The  wave  is  given.   The  lift  due  to  the  wave  can  be 
solved  explicitly  by  Section  B.   However,  the  required  con- 
trol surface  motion  (ALPHA  and/or  BETA)  is  only  implicitly 
expressed  in  Equations  (41)  and  (47)  of  Section  C.   Both 
Equations  (41)  and  (47)  are  complex  expressions.   Each  one 
represents  two  equations  either  as  the  real  and  imaginary 
parts  or  the  magnitude  and  phase  angle. 

The  strategy  used  to  solve  for  ALPHA  and/or  BETA  is: 

1,   Substitute  in  Equations  (41)  and/or  (47)  expressions 


for  ALPHA  and/or  BETA  and  their  derivatives  in  terms 
of  an  unknown  a0  and  <!>!.  (30/a0  and  A(J>  known  con- 
stants ) 

2.  Expand  multiple  angle  expressions. 

3.  Separate  the  real  and  imaginary  parts  of  the  modi- 
fied Equations  (41)  and/or  (47).  (-ReLIFTW, 
-ImLTFW) 

4.  Divide  out  cx0 ,  Cos<h,p,  and  b.  (-ReLIFTW/(a0Cos<{ilpb) . 

5.  Group  the  remaining  terms  into  terms  containing 
Tan  $i-- \,   Xlf  and  constant  terms — Y2>  X2. 
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-ReLIFTW/UoCos^pb)    =  Tan<J>i(Yi)   +    (Y2)  (50) 

-Im!TFW/(aoCos4»ipb)   =  Tan^iUJ   +   (X2)  (5D 

6.  Divide  the  real  equation  by  the  imaginary   equation. 

-RetTFW/(-ImIXFT¥)»   R 

-   (Tan4>l(Yl)+(Y2))/(Tan<|>l(Xl)   +   (X2))  (52) 

7.  Separate  out  Tan  <J>i 

Tan  0i  =  (Y2-RX2)/(RX,-Yl)  (53) 

8.  Solve  for  4>i. 

9.  Magnitude  a0  is 

a0  «  -ReLIFTW/[Cos<J>1pb(Tan(J>l(Y1)  +  (Y2))]     (5*0 

10.   The  problem  is  solved. 
See  Appendix  C  for  details.  With  ALPHA  and  BETA,  the  equa- 
tions in  Section  C  will  give  lift  and  moment  from  the  motion. 

E.   SUMMATION  OP  HYDRODYNAMIC  FORCES 

Summation  of  LIFTW,  LIFTA  and/or  LIFTS'  should  equal 
zero.   Summation  of  MOMTW,  MOMTA  and/or  MOMTB  give3  the 
total  hydrodynamic  moment  about  ba.   Summation  of  MOMTB W, 
M0MTBA  and/or  M0MTB3  gives  the  total  hydrodynamic  moment 
acting  on  the  flap  about  be. 
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F.   REQUIRED  APPLIED  MOMENTS 

The  control  system  must  apply  a  moment  (MOMAP)  on  the 
control  surface  to  produce  the  required  motion.   The  dyna- 
mics of  this  interaction  is  expressed  in  Equations  (55) 
and  (56). 


I  -ALPHADDOT  -  MOMTW  +  MOMTA  +  MOMTB  +  MOMAPA 


(55) 


Ib-BETADDOT  =  MOMTBW  -f  MOMTBA  +  MOMTBB  +  MOMAPB    (56) 

Prom  these  equations,  the  required  applied  moments  can  be 
found.   Figure  (21)  shows  these  relationships  in  phasor  form. 
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G.   AVERAGE  POWER 

The  average  power  (PAVE)  required  to  produce  the 
motion  is  expressed  in  Equations  (57)  and  (58). 

PAVEA  -  |  aTPHADOT| • | MOMAPA | -POWER  FACTOR/2         (57) 

PAVEB  =|BETADOT|- |MOMAPB | • POWER  FACTOR/2  (58) 

where  the  power  factor  is  the  cosine  of  the  phase  angle 
between  the  motion  velocity  and  the  applied  moment. 

H.   PROGRAM 

Two  programs  were  written  in  Fortran  IV  for  use  at 
the  Information  Processing  Center  on  the  "IBM  System  360". 
The  first  program  writes  out  the  details  of  the  power 
calculation,  and  the  second  writes  out  a  minimum  of  identi- 
fying Information  and  plots  the  power  results  on  a  "Calcomp 
Plotter".   An  outline  of  the  program  is  shown  on  Figure  (22). 

The  encounter  frequency  is  calculated  using  Equation 
(27)  as  shown  in  Figure  (10).   The  reduced  frequency  Is 
calculated  using  Equation  (20)  but  taking  the  absolute 
magnitude  for  use  with  Bessel  Functions.   The  Theodorsen 
Function  is  calculated  using  Equation  (31)  by  calling  the 
"IBM  Scientific  Subroutine  Package"  routines  BESJ  and 
BESY.   Because  of  the  iterative  method  of  calculation  used 
in  BESY,  the  reduced  frequency  parameter  k  can  not  be 
less  than  0.05.   This  is  guarded  against  in  the  program 
but  results  In  errors  in  the  Theodorsen  Function  for  reduced 
frequency  below  0.05. 
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The  WAVLIF  routine  calculates  the  complex  phasor 
forces  generated  by  the  wave  using  the  analysis  of  Section 
B.   The  lift  due  to  waves  (LIFTW)  is  brought  into  the  CONTRL 


routine.   Here  the  complex  phasor  ALPHA  and/or  BETA  are 
calculated  to  produce  an  equal  and  opposite  lift  to  cancel 
LIFTW.  The  procedure  used  is  described  in  Section  D.   The 


motion  (ALPHA  and/or  BETA)  is  brought  into  the  two  routines 
LIFA  and  LIFB  to  calculate  the  lift  caused  by  the  motion,  1) 
as  a  check  on  the  results  of  CONTRL;  and  2)  in  order  to 
calculate  the  moments  about  ba  and  be,  using  the  equations 
of  Section  C. 

So  far,  the  control  surface  motion  and  the  hydrodynamic 
moments  are  known.   The  hydrodynamic  moments  due  to  waves 
and  control  surface  motion  are  summed  as  in  Section  E.   Next, 
the  required  applied  moments  are  calculated  as  in  Section  F. 
The  average  power  is  calculated  using  the  analysis  of  Section 
G.   Note  all  results  are  2-D  per  foot  of  span.   An  important 
requirement  is  the  compatibility  between  mode  specification 
and  the  input  data  values  for  BA  ■  0o/ao,  and  DPHE  *  A<J>. 


For  Mode  1  ALPHA  motion, 

BA  »  0.0  DPHE  -  0.0 
For  Mode  2  BETA  motion, 

BA  »  1.0  DPHE  =0.0 
For  Mode  3  ALPHA  and  BETA  motion 

DPHE  *  0.0 
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The  plotting  routines  plot  the  motion  variable  magni- 
tude and  phase  angle,  and  the  average  power  as  functions 
of  encounter  frequency.   One  plot  per  foil,  the  plotter 
treats  each  wave  run  as  a  point  and  connects  each  point 
with  a  straight  line.   The  choice  of  wave  Input  data  should 
be  arranged  to  give  a  somewhat  continuous  variation  a3  was 
done  for  the  series  discussed  In  Section  I. 

Appendix  D  contains  a  listing  of  the  programs  and  an 
example  of  their  written  output.   Examples  of  the  plotted 
output  can  be  found  in  Section  I. 
I.   RESULTS 

The  two  pieces  of  input  data  required  for  each  run  of 
the  program  are  the  foil  data  and  the  wave  data.  Three  sets 
of  calculations  were  performed  for  the  three  candidate  foils, 
as  described  by  Table  (2): 

Mode  1  (Full  Incidence) 

Mode  2  (Trailing  Edge  Flap) 

Mode  3  (Tab  Foil) 

The  first  set  is  a  series  of  waves  with  ship  speed  con- 
stant, wave  length  constant,  and  angle  of  attack  (SI)  between 
the  ship  and  the  waves  varying  from  "dead  ahead"  to  "dead 
astern".   Table  (6)  displays  the  wave  parameters  for  each 
run  of  thi3  set.  Figures  (23),  (24),  and  (25)  are  the  compu- 
ter plotted  results  for  the  three  foils  run  with  this  series. 
The  series  can  be  pictured  as  a  slow  180°  turn  being  executed 
by  the  ship  at  constant  speed  in  regular  waves  of  constant 
wave  length. 
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TABLE  (6) 


SHIP  ANGLE  OF  ATTACK  TO  THE  WAVE,  VARIATION 


HW-n 

WL=X 

u=u 

SI** 

PHEW=<j>w 

RUN   0 

(FT) 

(FT) 

(Fr/SEC) 

(RADS) 

( RADS ) 

1 

7.0 

280,0 

80.0 

0.0 

-1.5708 

2 

7.0 

280.0 

80.0 

.3927 

-1.5708 

3 

7.0 

280.0 

80.0 

.7854 

-1.5708 

4 

7.0 

280.0 

80.0 

1.1781 

-1.5708 

5 

7.0 

280.0 

80.0 

1.3745 

-1.5708 

6 

7.0 

280.0 

80.0 

1.5708 

-1.5708 

7 

7.0 

280.0 

80.0 

1.7672 

-1.5708 

8 

7.0 

280.0 

80.0 

1.9635 

-1.5708 

9 

7.0 

280.0 

80.0 

2.1599 

-1.5708 

10 

7.0 

280.0 

80.0 

2.3562 

-1.5708 

11 

7.0 

280.0 

80.0 

2.5526 

-1.5708 

12 

7.0 

280.0 

80.0 

2.7489 

-1.5708 

13 

7.0 

280.0 

80.0 

2.9453 

-1.5708 

14 

7.0 

280.0 

80.0 

3.1416 

-1.5708 
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The  second  set  Is  a  series  of  waves  with  ship  speed  (U) 
varying  from  90  ft /sec  down  to  25  ft/sec,  wave  length  con- 
stant, and  angle  of  attack  between  the  ship  and  the  waves  at 
"dead  ahead"  and  "dead  astern".   Table  (7)  displays  the  wave 
parameters  for  each  run  of  this  set.  Figures  (26),  (27)* 
and  (28)  are  the  computer  plotted  results  for  the  three  foils 
run  with  this  series.   The  series  can  be  pictured  as  a  slow 
deceleration  by  the  ship  in  ahead  regular  waves  of  constant 
wave  length  and  a  slow  acceleration  by  the  ship  in  astern 
regular  waves  of  constant  wave  length. 

The  third  set  is  a  series  of  wave 3  with  constant  ship 
speed,  wave  length  (WL)  varying  from  280  ft  to  140  ft,  and 
angle  of  attack  between  the  ship  and  the  waves  at  "dead 
ahead"  and  "dead  astern".   Table  (8)  displays  the  wave  para- 
meters for  each  run  of  this  set.   Figures  (29),  (30),  and 
(3D  are  the  computer  plotted  results  for  the  three  foils 
run  with  this  series.   The  series  can  be  pictured  as  a  slow 
transit  by  the  ship  at  constant  speed  in  ahead  regular  waves 
of  reducing  wave  length  and  a  slow  transit  in  astern  regular 
waves  of  increasing  wave  length. 

Additional  trial  and  error  runs  were  made  in  order  to 
discover  a  useful  motion  ratio  30/a0  (BA)  and  a  useful  phase 
angle  difference  A<j>  (DPHE)  for  Mode  3  (Tab  Foil).   The  search 
is  not  complete,  but  BA  ■  2.0  and  DPHE  -  -.3927  radians 
appeared  best  and  were  utilized  for  the  Mode  3  runs. 
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TABLE  (7) 


SPEED  VARIATION 


hw»ti 

WL-X 

U«U 

31=^ 

PHEW=4»w 

RUN    if 

(FT} 

(FT) 

(FT/SEC) 

(RADS) 

(RADS) 

1 

7.0 

280.0 

90.0 

0.0 

-1.5708 

2 

7.0 

280.0 

85.0 

0.0 

-1.5708 

3 

7.0 

280.0 

80.0 

0.0 

-1.5708 

4 

7.0 

280,0 

75.0 

0.0 

-1.5708 

5 

7.0 

280.0 

70.0 

0.0 

-1.5708 

6 

7.0 

280.0 

65.0 

0.0 

-1.5708 

7 

7.0 

280.0 

60.0 

0.0 

-1.5708 

8 

7.0 

280.0 

55.0 

0.0 

-1.5708 

9 

7.0 

280.0 

50.0 

0.0 

-1.5708 

10 

7.0 

280.0 

45.0 

0.0 

-1.5708 

11 

7.0 

280.0 

40.0 

0.0 

-1.5708 

12 

7.0 

280.0 

35.0 

0.0 

-1.5708 

13 

7.0 

280.0 

30.0 

0.0 

-1.5708 

m 

7.0 

280.0 

25.0 

0.0 

-1.5708 

15 

7.0 

280.0 

25.0 

3.1416 

-1.5708 

16 

7.0 

280.0 

30.0 

3.1416 

-1.5708 

17 

7.0 

280.0 

35.0 

3.1416 

-1.5708 

18 

7.0 

280.0 

40.0 

3.1416 

-1.5708 

19 

7.0 

280.0 

45.0 

3.1416 

-1.5708 

20 

7.0 

280.0 

50.0 

3.1416 

-1.5708 

21 

7.0 

280.0 

55.0 

3.1416 

-1.5708 

22 

7.0 

280.0 

60.0 

3.1416 

-1.5708 

23 

7.0 

280.0 

65.0 

3.1416 

-1.5708 

24 

7.0 

280.0 

70.0 

3.1416 

-1.5708 

25 

7.0 

280.0 

75.0 

3.1416 

-1.5708 

26 

7.0 

280.0 

80.0 

3.1416 

-1.5708 

27 

7.0 

280.0 

85.0 

3.1416 

-1.5708 

28 

7.0 

280.0 

90.0 

3.1416 

-1.5708 
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TABLE  (8) 


WAVE  LENGTH  VARIATION 


hw-ti 

WL»X 

U=U 

SI*ij; 

PHEW*<J>w 

RUN   it 

(FT) 

(FT) 

(FT/SEC) 

(RADS) 

( RADS ) 

1 

3.5 

140.0 

80.0 

0.0 

-1.5708 

2 

4.2 

168.0 

80.0 

0.0 

-1.5708 

3 

M 

196.0 

80.0 

0.0 

-1.5708 

4 

5.6 

224.0 

80.0 

0.0 

-1.5708 

5 

6.3 

252.0 

80.0 

0.0 

-1.5708 

6 

7.0 

280.0 

80.0 

0.0 

-1.5708 

7 

7.0 

280.0 

80.0 

3.1416 

-1.5708 

8 

6.3 

252.0 

80.0 

3.1416 

-1.5708 

9 

5.6 

224.0 

80.0 

3.1416 

-1.5708 

10 

*».9 

196.0 

80.0 

3.1416 

-1.5708 

11 

4.2 

168.0 

80.0 

3.1416 

-1.5708 

12 

3.5 

140.0 

80.0 

3.1416 

-1.5708 

55 


ex 


cr 


<x 

X 
0_ 


cc 

cr 

X 


v. 

03 


a 


cc 


— I — 


■  1  .bCi 


UT 


if  C1U 


"I — 


56 


t. 


D 

-ZL 

figure    23 

Full 

Incidence 
(SI)  Var. 

A=  la| 
Z  =  La 
P  =  pavea 


57 


v 


■Y 


CD 


aa  "7 


03   T 


an 


C3 


ar 


a: 


— I — 


kE  m 


— r 


'0 


— r~ 


58 


' 


t) 


D 


CD 


figure  24 

Trailing 
Edge  Flap 
(SI)  Var 

B  =  \p\ 
Y  =LP 

Q  =  paveb 


59 


2L 
CL 


CD 
CL 


CO      ' 


C3 


az 


i 


5  <?b 


60 


— I — 


CD 


in 


_Q 


v~> 


CD 


figure  25 

Tab 

Foil 

(SI)  Var. 

A  = 

Ia| 

\p\  = 

2-lal 

Z  = 

/a 

LP  = 

La-  22* 

P  = 

pavea 

Q  = 

- paveb 

61 


CD 


cr 


cr 

cr 

— 

cr 


CD 


a 


i 


-2  .25 


il 


62 


^z 


-^. 


figure  26 
Full 

Incidence 
(U)  Van 

A=  Ia| 
Z  =  La 
P  =  pavea 


63 


or 


CD 


Ll_ 


3 


61 


OQ 


figure  27 
Trailing 
Edge  Flap 
(U)  Var 

B  =  l/?l 

Y=/_£ 

Q  =  paveb 


65 


11 


-z 


CD 

cr 


2: 
cr 

cr 

3Z 
Q_ 

I 

cr 


CD     ' 


C3 


2 


1 


-2   2b 


a 
ex 


.25 


66 


21 


x. 


m 


figure  28 

Tab 
Foil 
(U)  Van 

A=la| 
l/?l  =  2lal 

Z  -1st. 
L/3  =/a  -22'/z 
P  =  pavea 
Q  = - paveb 


67 


X 


CD 


CL 


C3 

cr 

x 

CL 


CD 
CD 


a 

a 


-ij  ^ 


i — 


a 
cr 


68 


21 


*S\ 


figure  29 
Full 

Incidence 
(WL)  Van 

A  =  la| 
Z  =  la. 
P  =  pavea 


69 


T 


tn  7 


z: 


en 

CD 


cc 


! 


"""1 


-  .75 


kE   RRD/^EC 


^ 


CD 


~r- 


5    75 


1 — 


70 


O 


^ 


■ 


figure  30 
Trailing 
Edge  Flap 
(WL)Var 

B=  \p\ 
Y  =  LP 
Q  -  paveb 


71 


CD 


Q_ 


CD 
T. 


1 


CD 


C3 
CL. 


cc  S 


CL, 


tr 


RRD/ 


7 


72 


lO 


o 

CD 


CD 
CD 


CD 


ZD 


- 


10 

figure  31 
Tab 
Foil 
(WL)  Var. 

A  =  lal 
l/?l=  2-lal 

Z   =/a 
Z./?  =/A  -22^ 

P  =  pavea 
Q  =  -  paveb 


73 


All  the  results  are  in  error  for  the  range  of  reduced 
frequency,  -,0  5<k<+.05  »  because  of  the  fix-up  prooedure 
mentioned  in  Section  H.   The  error  is  dramatic  because  of 
the  sensitivity  of  power  to  the  related  phase  angles  of  its 
components.   The  error  is  zero  at  k  *  0,  and  beyond  the 
range  k  =±.05.   The  range  for  k  translates  into  a  differ^ 
ent  range  of  u>e  for  each  of  the  three  wave  series.   The 
reason  is  that  k  is  composed  of  both  u>  and  U,  which  are 
not  solely  functions  of  one  another,  Equation  (27).   For 
a  given  U,  u>e  remains  dependent  on  X  and  ^.   The  following 
are  the  ranges  of  u>e  which  are  in  error  for: 

Series  1.   (SI)  Variation 

-.85<w  <+.85    (Rads/Sec) 
Series  2.   (U)  Variation 

-.94<we<+.29    (Rads/Sec) 
Series  3.   (WL)  Variation 

Runs  not  in  the  range 

The  plotted  results  are  somewhat  misleading  in  a  part 
of  the  u>e  range,  because  the  plotter  plots  straight  lines 
between  points.   The  following  are  the  ranges  in  a>e  where 
there  are  significant  gaps  in  plotted  points: 

1.   (SI)  Variation 

None 
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2.   (U)  Variation 


+  .28<u>e<+l.Ml    (Rads/Sec) 

3.   (WL)  Variation 

-,9^<w  <+2.6*i    (Rads/Sec) 
e 

Aside  from  the  fix-up  error  and  taking  note  of  the  range  for 
the  gap  in  points,  the  plotted  results  give  a  reasonably 
good  picture  of  the  power  required. 

In  the  case  of  Mode  3  ALPHA  and  B"ET"A  Motion  (Tab  Foil), 
the  plotter  plots  the  Beta  power  as  Q=»  -  PAVED.   This  opera- 
tion places  the  Q  result  relative  to  the  P  result  (PAVEA) 
such  that  the  distance  between  them  is  the  total  power  re- 
quired by  the  system.   If  the  distance  between  Q  +  P  is  in 
the  positive  direction,  the  total  power  required  is  positive. 
If  the  distance  between  Q  +  P  is  in  the  negative  direction, 
the  total  power  required  is  negative .   The  above  is  true 
regardless  of  where  the  plots  lie  on  the  absolute  scale. 

THE  FOLLOWING  ARE  COMMENTS  ON  EACH  FOIL  RESULT: 

1.  Alpha  Motion  (Full  Incidence)  Mode  1 

a.   (SI)  Variation*     Figure  (23)  shows  a  decreasing 
average  power  requirement  from  "dead  ahead"  around  to 
quartering  seas  where  the  encounter  frequency  goes  to  zero 
(power  goes  to  zero);  and  then  continuing  around  to  "dead 
astern",  the  power  increases  again  for  we  ■  (-).   The  maximum 
power  occurs  at  SI  ■  0  and  it  is  positive. 
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The  magnitude  of  the  Alpha  motion  is  constant,  because 
the  disturbance  is  unchanged  in  magnitude  and  the  ship's 
(lift  generating)  speed  is  constant.   Note,  at  «  ■  0 
(SI  «  112°),  the  magnitude  of  Alpha  does  not  go  to  zero. 
This  is  because,  depending  upon  what  the  phase  angle  of  the 
wave  is  that  the  ship  is  running  along  (crest,  back- face, 
trough,  or  front-face),  Alpha  will  have  to  be  a  constant  0, 
(  +  >.  0,  (-). 

The  phase  angle  of  Alpha  relative  to  the  wave  height 
changes  from  lagging  90°  +  e  for  to  ■  (+),  to  leading  90°  -  e 
for  we  ■  (-).   At  u)e  ■  0,  the  phase  angle  is  90°  from  the 
wave  height;  but,  more  important,  it  is  180°  out  of  phase 
with  the  wave  upwash  velocity  UPW  as  is  expected  from  quasi- 
steady  analysis.   The  basic  characteristic  of  shifting  from 
lagging  to  leading  and  being  out  of  phase  with  UPW  at  we  ■  0 
is  common  to  all  the  foils  in  all  the  wave  series. 

b.   (U)  Variation.     Figure  (26)  show3  increasing 
magnitude  for  the  required  Alpha  motion,  as  speed  drops  from 
90  ft/sec  In  ahead  seas  (lesser  w  a  (  +  )).  The  magnitude 
decreases  as  the  ship  speed  increases  In  astern  seas  (greater 

we  -  <-». 

The  phase  angle  goes  through  a  similar  variation  as  in 
the  SI  series. 

The  average  power  decreases  as  the  ship  slows  in  ahead 
seas.   In  astern  seas,  the  power  dips  down  through  zero  at 
we  *  0  and  then  climbs  back  up  as  speed  Increases.   The  gap 
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between  points  at  ±  25  ft/sec  is  necessary  because  of  the 
limited  capability  of  the  lifting  foils  to  support  the 
ship  at  low  speed.   This  is  evident  by  the  rapidly  increasing 
motion  magnitude  (shown  in  Figure  (26))  required  to  produce 
the  lift  at  low  speeds. 

c.   (WL)  Variation.   Figure  (29)  shows  the  magnitude 
of  Alpha  motion  increases  slightly  as  the  wave  length 
increases  in  ahead  seas.   This  is  because  the  increased 
wave  length  presents  an  increasing  disturbance.  The  same 
variation  appears  in  astern  seas. 

The  phase  angle  changes  in  a  similar  fashion  as  in  the 
SI  series. 

The  average  power  decreases  as  wave  length  Increases 
in  ahead  seas  (lesser  u>e  «=  (  +  )).   In  astern  seas,  the  power 
again  decreases  as  wave  length  increases  (lesser  u>e  *  (-)). 
Note,  there  is  no  effect  of  C(k)  error  in  this  plot,  because 
u>e  is  out  of  the  range. 

2.  Beta  Motion  (Trailing  Edge  Flap)  Mode  2 

a.   (SI)  Variation.     Figure  (24)  shows  the  magnitude 
of  the  Beta  motion  constant  as  in  the  Alpha  case.   The  magni- 
tude of  Beta  is  larger  than  the  Alpha  case  which  is  to  be 
expected  for  a  1/k   chord  flap  compared  to  full  incidence. 

The  phase  angle  goes  through  a  similar  variation  as  in 
the  Alpha  case. 
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b.  (U)  Variation.     Figure  (27)  shows  the  magnitude 
increasing  in  a  similar  manner  as  in  the  Alpha  case.   Again, 
the  magnitude  of  Beta  is  larger  than  that  for  Alpha. 

The  phase  angle  change  is  similar  to  the  (SI)  variation. 

The  average  power  follows  a  similar  trend  as  in  the 
Alpha  case;  but,  again,  it's  sensitivity  to  error  in  C(k) , 
for  the  range  -.94<w  <+.29,  is  much  greater.   The  only  point 
that  is  correct  in  this  range  is  at  u>e  ■  0.   Again,  it  is 
speculated  that,  if  correct  values  for  C(k)  were  used,  the 
power  curve  would  fair  into  zero  at  u  ■  0  and  climb  back 
up  again  in  a  similar  fashion  to  the  Alpha  case. 

c.  (WL)  Variation.     Figure  (30)  shows  the  magnitude 
of  Beta  increasing  in  a  similar  fashion  as  in  the  Alpha  case. 

The  phase  angle  changes  in  the  same  manner  as  in  the 
Alpha  case. 

The  average  power  increases  for  increasing  wave  length 
in  ahead  seas  (lesser  we  *  (  +  )).   This  result  is  in  contrast 
with  the  Alpha  result  where  power  decreased  for  increasing 
wave  length  in  ahead  seas.   The  average  power  in  astern  seas 
decreases  for  increasing  wave  length  (lesser  we  *=  (-)). 
Note,  there  is  no  effect  of  C(k)  error  in  this  plot,  because 
«e  is  out  of  the  range. 

3.   Alpha  and  Beta  Motion  (Tab  Foil)  Mode  3 

a.   (SI)  Variation.     Figure  (25)  shows  the  magnitude 
of  Alpha  as  a  constant  straight  line  as  in  the  Alpha  only 
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case.  The  associated  Beta  magnitude  is  two  tines  the  Alpha 
magnitude;  however,  it  is  not  plotted.   This  is  true  for  the 
(U)  variation  and  (WL)  variation  as  well. 

The  Alpha  phase  angle  has  a  similar  trend,  as  before, 
but  the  o)e  »  (±)  variations  are  displaced  away  from  90°. 
The  Alpha  motion  now  lags  the  wave  height  by  90°  -  e  for 
we  ■  (  +  ),  and  leads  by  90°  +  e  for  w  »  (-).  The  ramp-type 
Jog  near  u>e  ■  0  is  caused  by  the  plotter  moving  from  we  = 
+.16  ■*  -.14.  The  Jog  would  approach  a  step  type  for  points 
approaching  «e  ■  0.   The  Beta  phase  angle  lags  the  Alpha 
phase  angle  by  22  1/2°  for  both  u>  ■  (±);  however,  it  is  not 
plotted.  This  is  true  for  the  (U)  variation  and  (WL)  varia- 
tion as  well. 

The  results  for  power  are  important.   As  was  pointed  out 
earlier  in  this  section;  if  the  distance  between  Q  (Beta 
power)  ♦  P  (Alpha  power)  is  in  the  negative  direction,  the 
total  power  to  drive  the  motion  is  negative  (the  control 
system  must  absorb  power  from  the  motion).   Figure  (25)' s 
plot  of  Q,  P  is  exactly  this  case.   Figure  (32)  demonstrates 
in  qualitative  form  the  relative  position  of  the  motion 
phasors  and  their  lift  phasors. 
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b.  (U)  Variation.     Figure  (28)  shows  the  magnitude 
of  Alpha  increasing  as  in  the  Alpha  only  case. 

The  Alpha  phase  angle  variation  is  similar  to  the 
Mode  3  SI  case. 

The  average  total  power  is  negative  over  the  range  of 
u)e  except  near  u)e  ■  +1.50.   Here  the  power  required  by  Alpha 
is  positive,  but  the  power  required  by  Beta  is  Just  enough 
negative  for  the  total  power  to  be  zero. 

c.  (WL)  Variation.     Figure  (31)  shows  a  slightly 
increasing  Alpha  magnitude  as  in  the  Alpha  only  case. 

The  Alpha  phase  angle  is  changing  in  a  manner  similar  to 
the  SI  case. 
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The  total  average  power  starts  out  positive  for  small 
wave  lengths  in  ahead  seas.   At  u>e  »  +4.28,  the  total  power 
is  zero  and  then  goes  negative  for  decreasing  coe.   The 
average  negative  power  for  Beta  is  not  enough  negative  to 
cancel  the  positive  Alpha  power  for  we  +4.28.  Note,  the 
C(k)  error  does  not  effect  this  plot. 

From  the  above  results,  it  is  clear  that  an  under- 
standing of  the  unsteady  effects  on  foils  is  not  only 
necessary  but  also  rewarding.   Table  (9)  lists  some  of  the 
important  power  results. 


SI 


U 


WL 


TABLE  (9) 

(Full  Incidence)  (Trailing  Edge  Flap)  (Tab  Foil) 

Alpha  Motion       Beta  Motion  Alpha  and  Beta  Motion 

Maximum  ±         Maximum  ±  Maximum  ± 

Average  Power      Average  Power  Average  Power 

ft-lb/sec          ft-lb/sec  ft-lb/sec 

per  ft  span per  ft  span  per  ft  span 

+  738.5          +  468.0  -  15^.9 

+  774.5                         +  534.7  -  153.5 

+1105.2           +H91.2  +  23.6,  -  133.8 
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V.   DRAG 

The  components  of  drag  for  a  foil  are: 

Skin  Friction 

Pressure  Drag 

Induced  Drag 

Wave  Drag 

Interference  Drag 
Each  of  these  has  a  contribution  from: 

Steady  Motion 

Waves 

Foil  Motion 
The  task  of  sorting  out  all  of  these  parts  is  formidable. 
For  a  2-D  analysis,  the  induced  parts  can  be  eliminated. 
For  a  small  strut,  the  interference  drag  can  be  neglected. 
For  a  well-submerged  foil,  the  wave  drag  can  be  neglected. 
Assuming  small  motions  with  no  separation,  cavitation  or 
ventilation,  the  pressure  drag  can  be  neglected. 

This  leaves  the  effects  of  steady  motion,  waves,  and 
foil  motion  on  skin-friction  drag.  The  temptation  is  to 
linearly  superpose  these  three  effects  as  was  done  for  lift. 
Yet  the  model  for  skin  friction  is  the  boundary  layer. 

Reference  (2)  describes  the  effect  on  foil-skin  friction 
by  a  uniform  inflow  with  a  harmonic  perturbation  superim- 
posed.  This  is  a  model  of  the  effect  of  waves  on  the  skin 
friction  of  a  steady  foil.   The  shearing  stress  lead3  the 
velocity  perturbation  (horizontal  orbital  velocity  for  the 
wave),  and  the  phase  angle  approaches  4 5°  for  k  -*•  ». 
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An  oscillating  cylinder  in  a  fluid  at  rest  is  discussed, 
and  an  oscillating  infinite  flat  plate  in  a  fluid  at  rest  is 
discussed.  The  problem  of  a  finite  chord  oscillating  flat- 
plate  foil  in  uniform  flow  is  not  discussed.   No  mention  is 
made  of  the  advisability  of  superposing  such  a  solution  on 
the  effects  due  to  waves. 

Although  the  required  motion  of  the  control  surface 
results  in  zero  total  lift,  the  lift  distribution  is  differ- 
ent for  the  wave  and  for  the  motion.  This  is  evidenced  by 
the  different  values  for  wave  moments  and  motion  moments. 
Therefore,  the  pressure  distribution  will  oscillate  back 
and  forth  across  the  foil.  It  is  speculated  that  the  problem 
can  no  longer  be  solved  by  linear  superposition,  even  if 
there  v;as  a  solution  to  the  foil  motion  effects.   The  solu- 
tion of  the  coupled  effects  of  waves  and  motion  appears 
beyond  analytical  solution  today. 

All  that  is  left  is  skin-friction  drag  from  steady 
motion  in  uniform  flow.  The  result  for  the  three  candidates, 
with  the  same  overall  foil  geometry,  is  that  there  Is  no 
significant  difference  in  skin-friction  drag. 

This  statement  is  adequate  for  comparing  the  overall 
power  of  the  three  candidates.  However,  if  comparisons  of 
these  foils  were  being  made  with  other  systems,  such  as 
spoilers,  the  above  rationalization  is  not  adequate  because 
ventilation  and  separation  play  a  much  different  part  in  the 
drag  of  such  systems. 
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VI,   WEIGHT 

Any  statement  about  weight  without  a  first  iteration 
of  the  complete  control  system  design  must  be  qualitative. 
However,  for  the  effect  of  weight  to  influence  the  choice 
of  control  system,  an  absolute  figure  must  be  found.  For 
the  purpose  of  emphasizing  that  power  must  be  supplied  to 
lift  the  weight  of  the  control  system,  an  approximate 
analysis  will  be  pursued. 

A.  REFERENCE  FULL-INCIDENCE  SYSTEM 

To  use  as  a  reference  system,  the  weight  of  a  full- 
incidence  (Alpha  motion)  control  system  is  estimated  as 
follows: 

Group  I    Foil  Weight  »15,700  lbs  *  7*»0  lbs/ft 

Group  II  Mechanical  Drive  Links     «   650  lbs  *  30  lbs/ft 

Hydraulic  Servo  Cylinders  «   650  lbs  *  30  lbs/ft 

Hydraulic  Fluids  and  Pipes  ■  1,070  lbs  «  50  lbs/ft 

Hydraulic  Pumps  and  Drive  «   ^00  lbs  «  20  lbs/ft 

Group  III  Autopilot  and  Computer     «   200  lbs  *  10  lbs/ft 

Group  IV  Sensors  ■   200  lbs  «   10  lbs/ft 

B.  COMPARATIVE  WEIGHT 

The  Group  I,  III,  and  IV  weights  for  the  three  candidates 
will  not  be  significantly  different.   The  Group  II  weights 
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will  approximately  vary  linearly  with  average  power.   For 
comparison,  from  Table  (9): 

Trailing  Edge  Flap  Max  Power   B   53^.7       Uo 
Full  Incidence  Max  Power        1105.2      *M0 


Tab  Foil  Max  Power        »  ~15|[.9       ,j. 
Full  Incidence  Max  Power        1105.2 

(Here  the  maximum  negative  power  is  taken  to  size  the 
servo  machinery  even  though  it  must  absorb  the  power. ) 


C.  TRAILING  EDGE  FLAP  SYSTEM  WEIGHT 

Group  II   *  70  lbs/ft 

D.  TAB  FOIL  SYSTEM  WEIGHT 

Group  II   *  20  lbs/ft 

E.  POWER 

The  weight  becomes  power  through  the  lift/drag  ratio 
and  the  speed  of  the  hydrofoil  ship: 

Power  =  (Drag/Lift) 'Weight 'Speed  (59) 

For  comparative  purposes,  the  lift/drag  ratio  is  7  and  the 
speed  is  80  ft/sec.   Power  results  due  to  weight  are  as 
follows: 

Full  Incidence  =*  1500  ft-lb/sec  per  ft  span 
Trailing  Edge  Flap  «  800  ft-lb/sec  per  ft  span 
Tab  Foil  «   230  ft-lb/sec  per  ft  span 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

In  summary,  the  contributions  of  power,  power  due  to 
drag,  and  power  due  to  weight  have  been  developed  for  three 
candidate  hydrodynamic  devices  for  use  in  a  constant  lift, 
wave  alleviation,  control  system  for  a  hydrofoil  ship. 

A.   OVERALL  CONTROL  SYSTEM  POWER 

A  3ervo  system  efficiency  is  required  to  convert  the 
control  surface  power  into  input  power  to  the  system.   An 
efficiency  of  .6  is  assumed  and  has  been  applied  to  the 
maximum  power  results  in  Table  (9)» 

Table  (10)  shows  the  results  for  overall  control  system 
power. 

TABLE  (10) 


Full  Incidence 
ft-lb/sec 
per  ft  span 


Trailing 
Edge  Flap 
ft-lb/sec 
per  ft  span 


Tab  Foil 
ft-lb/sec 
per  ft  span 


Power 


+  18U0. 


+  890. 


+  40. 


Drag 
Power 


Weight 
Power 


+  1500. 


800. 


+  230. 


Total 


+  33*»0. 


+  1690. 


+  270. 
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From  this  analysis,  the  tab  foil  looks  very  promising  for 
use  as  the  hydrodynamic  device  in  a  constant  lift  control 
system.   Development  of  the  sensor,  servo,  and  autopilot 
system  designs  for  use  with  the  tab  foil  is  recommended. 

If  the  choice  is  restricted  to  either  full  incidence 
or  trailing  edge  flap,  the  trailing  edge  flap  looks  most 
promising.  The  development  of  a  control  system  utilizing 
a  full  incidence  foil  for  stability  and  maneuverability , 
and  a  trailing  edge  flap  for  wave  alleviation,  is  recommend- 
ed. 

But  the  most  important  recommendations  are: 

1.  Wave  alleviation  should  be  an  important  considera- 
tion in  the  design  of  an  open-ocean  hydrofoil  ship. 

2.  Unsteady  lift  analysis  should  be  used  during  the 
design  of  the  control  system  for  an  open-ocean 
hydrofoil  ship. 

B.  SENSORS 

The  control  system  must  be  able  to  sense  the  magnitude 
of  the  upwash  velocity,  the  direction,  and  the  encounter 
frequency  of  the  wave  the  foil  is  about  to  see*  enough  in 
advance,  in  order  to  control  proper  foil  motion.   Wave  height 
information  is  not  enough.   It  is  suggested  that  wave  height 
and  orbital  velocity  be  measured  ahead  and  behind  the  ship. 

C.  MOTION 

Figure  (35)  shows  full  incidence,  trailing  edge  flap, 
and  tab  foil  motion.   The  amplitudes  have  been  exaggerated, 
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but  their  relative  magnitudes  have  been  maintained.   The 
most  interesting  motion  is  that  of  the  tab  foil.   With 
22  1/2°  lagging  phase  angle  for  Beta  (-A<fr),  the  Alpha  motion 
i3  caused  by  the  moments  generated  by  the  Beta  motion.   The 
motion  appears  "natural"  in  contrast  to  the  22  1/2°  leading 
system  shown  below  it.   Figure  (36)  shows  motions  for  90° 
lagging,  180°  out  of  phase,  and  90°  leading. 
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Of  these,  only  the  90°  lagging  appears  "natural".   Preliminary 
checks  with  the  computer  program  for  80  ft/sec  speed  in 
limiting  ahead  seas,  indicate: 

90°  Lagging  }   Alpha  Power  <<  0 
Beta  Power  <<  0 
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180°  Out  of  Phase  >  0  <<  Alpha  Power 

Beta  Power  <  0 

90°  Leading   }  0  <<  Alpha  Power 

0  <<  Beta  Power 

The  "natural"  looking  motions  require  negative  power.   The 
other  motions  require  positive  power.  References  (3)  and 
(4)  discuss  the  use  of  a  tab  to  initiate  aileron  motion.   The 
tab  is  moved  in  the  opposite  direction  from  the  direction  the 
aileron  will  move.   At  first,  this  sounds  like  the  aileron 
tab  3hould  be  180°  out  of  phase.   However,  the  Beta  motion 
must  lead  (in  time)  to  cause  the  motion.   This  lead  (in  time) 
in  the  opposite  direction  is  a  lag  relative  to  the  aileron 
motion. 

For  the  hydrofoil  ship,  lags  near  90°  result  in  large 
negative  powers  to  be  absorbed  (like  the  flutter  situation). 
If  small  lags  around  22  1/2°  are  used,  small  negative  powers 
result. 

D.   EFFECTS  ON  HYDROFOIL  SHIP  OPERATION 

In  the  limiting  sea  (ahead  or  astern),  slowing  down  the 
ship  will  reduce  the  control  system  power  but  will  increase 
the  magnitude  of  the  control  surface  motion.   Since  most 
foil  arrangements  have  a  mechanical  limit  to  foil  travel, 
slowing  down  may  not  be  the  best  operational  change.   Running 
with  the  wave  requires  no  foil  motion  and  no  control  system 
power.   Running  with  the  crest  or  trough  requires  no  control 
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surface  angle.   Running  with  the  face  or  back  requires  full 
-  or  +  control  surface  angle.   Running  with  the  crest  is 
recommended  because  no  control  surface  angle  is  required  and 
maximum  submergence  of  the  foils  is  obtained. 

The  angle  of  attack  to  the  sea  will  reduce  the  control 
system  power  but  not  the  foil  motion  magnitude.   When  running 
along  the  wave,  the  face  or  back  of  the  wave  present  no  span- 
wise  flow;  but  the  control  surface  will  be  at  full  -  or  + 
angle.  Running  along  the  crest  or  trough  will  present  a 
small  spanwlse  flow  to  the  foils  and  struts,  but  no  control 
surface  angle  will  be  required.   Running  along  the  trough  is 
recommended  because  no  control  surface  angle  is  required  and 
maximum  submergence  of  the  foils,  on  the  beam,  is  obtained. 

Finding  smaller  waves  decreases  control  surface  motion 
in  all  cases,  and  will  reduce  the  control  system  power  for 
a  trailing  edge  flap  system  in  ahead  seas.   But,  smaller 
waves  increase  control  system  power  for  the  full  incidence 
system  in  both  ahead  and  astern  seas.   Smaller  waves  increase 
control  system  power  for  the  trailing  edge  flap  system  in 
a3tern  seas. 

It  is  recommended  that  wave  alleviation  unsteady  effects 
be  analyzed  before  determining  operational  policy. 

E.   STATISTICAL  ANALYSIS 

The  foil  responds  to  wave  upwash  velocity  as  an  input. 
Therefore,  an  orbital  velocity  spectrum  with  direction  infor- 
mation is  required.   This  is  a  3-dimensional-surface  from 
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statistical  wave  data.   Correspondingly,  the  foil  response 
spectrum  is  no  longer  a  single  curve  but  is  a  3-cllinens ional- 
surface  for  a  given  ship's  speed.   The  complete  picture  is  a 
series  of  response  surfaces  for  a  range  of  speeds. 

It  is  recommended  that  an  attempt  be  made  at  bringing 
the  unsteady  foil  response  analysis  into  a  statistical  for- 
mat.  From  this,  a  complete  and  realistic  description  of 
hydrofoil  ship  energy  requirements  could  be  made. 
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APPENDIX   A:      HYDRODYNAMIC   DEVICES    FOR  SECTION  III-A 

I.  ANGLE   OF  ATTACK   CHANGERS 

A.   TRANSLATING  FOIL  (SUSPENDED  FOIL) 

The  foil  vertical  velocity  must  match  that  of  the  incom- 
ing wave.   If  the  speed  of  the  foil  is  zero,  w  ■  ww.   The 
vertical  velocity  of  the  incoming  wave  will  be  |n|  ■  nww. 
The  velocity  of  the  foil  must  be  |z|  ■  |n|;  but,  if  u>e  *  u>w, 
than  | as |  ■  |n|.   In  other  words,  the  foil  motion  would  exact- 
ly follow  the  wave  surface.   If  speed  is  increased,  the 
magnitude  of  we  increases  and  |z|  <  |n|.   Such  large  excur- 
sions would  probably  not  be  feasible.   A  patent  for  a  similar 
scheme  is  in  Reference  (5). 

II.  CAMBER  CHANGERS 

A.  PARTIAL  CHORD  FLUID  (JET  FLAP) 

Much  work  has  been  done  on  Jet  flaps  for  high  \lft  on 
aircraft.   The  feasibility  of  pumping  water  through  the  foil 
has  not  been  determined  as  yet.   Fluidic  controls  could  be 
employed  to  change  the  direction  of  the  jet  forfwave  allevia- 

tion.   See  References  (6)  and  (7). 

'■•%■ 

B.  UPPER  SURFACE  MECHANICAL  (MEMBRANE) 

This  scheme  would  employ  a  second  skin  for  the  .upper 
surface  of  the  foil.   The  space  between  could  be  pressurized 
to  change  the  membrane  shape.   The  upper  surface  lift  load 
would  have  to  be  taken  in  tension  in  the  membrane.   The 
structure  appears  to  be  the  greatest  problem. 
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C.   UPPER  SURFACE  FLUID  (FLUID  CAMBER) 

The  concept  of  using  a  layer  of  fluid  on  the  upper 
surface  to  change  the  apparent  thickness  to  the  oncoming 
flow  was  pursued  to  some  detail  early  in  the  project.   The 
following  is  an  approximate  analysis  for  such  a  scheme. 

In  order  to  change  the  lift  coefficient  of  the  foil 
from  .J*  to  .9,  the  camber  must  change  from  2.2#  of  chord  to 
5.2JC  of  chord,  or  a  change  of  3%   of  chord.   In  order  to 
change  the  camber  by  3%,   the  thickness  must  be  changed  by  6% 
of  chord. 

The  proposed  scheme  of  using  "fluid  camber"  to  change  a 
foil's  lift  characteristic  calls  for  pumping  fluid  from  the 
leading  edge  on  one  surface  of  the  foil  and  retrieving  the 
same  fluid  back  into  the  foil  at  the  trailing  edge. 
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After  some  fruitless  searching  for  a  potential  flow 
model  for  the  concentrated  Jet  at  the  leading  edge,  a  half- 
plane  source  and  sink  were  chosen  as  a  simple  model  of  both 
entry  and  exit  of  the  flow.   The  stream  function  for  a  source 
and  sink  of  equal  strength  m  and  separated  by  a  distance  d  is: 

-Ji  «  m  m^-i    -2y(l/2)d 
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The  stream  function  for  a  uniform  flow  in  the   x  direction 
at  speed  Uo  is: 

*  -  U0y  (2) 

For  the  superposition  of  these  flows: 

*  -  U.y  +  ^Tan-  ^ffjjffifa.  (3) 

We  need  the  half  plane  thickness   of  the  ellipse   formed  at 
ty  -  0  to  be  6%  of  chord,   or  y/d  ■   .06.      With  ij»  ■  0,   x  -  0, 
y  *    .  06d,   we   can  solve   for     m. 


TT7       Tan-iC?an  B)  (H) 


where 


Tan(B)   *       -2(.06)(l/2)d 
Tan{B)        (,06d)t-((i/2)d)» 


(5) 
Tan(B)   -   .243  (6) 

B  «  -166.67°  igff  (7) 

m  -    .1295«d.U0  (8) 

But  we  need  only   1/2  this   flow  but  across   the   span  of  the 
foil. 


Sp.  *   .0648*   tfo.ds  (9) 


Q  «    ,0648*U0»A  ft'/sec  (10) 
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For  an  ocean  hydrofoil  ship  of  300  tons  at  40  knots 
with  lift  coefficient  of  .*!,  requires  this  area: 


A   (I/2)pV2C; 


*  m   300*22^0      .  o68  ft* 

A      TT7TPCtrtW^WCJT      3b0  rt  (12) 

This  gives  us  a  required  flow  rate  for  the  "fluid  camber" 
system  of: 

Q  »  .0618  (67.5X368)  »  1610  ft'/sec 

or  =  7.22  x  10$  gal/min   (13) 

The  "fluid  camber"  scheme  is  not  feasible  for  the  full  scale 
hydrofoil  ship. 

III.   CAMBER  AND  ANGLE  OF  ATTACK  CHANGERS 

A.  ROTATING  FOIL/PARTIAL  CHORD  FLUID  (JET  TAB) 

This  scheme  would  utilize  the  jet  momentum  and  effect  on 
the  overall  flow  to  force  the  Alpha  motion.   Reference  (7) 
discusses  such  a  scheme.   Hopefully,  the  flow  rate  required 
in  the  Jet  could  be  reduced  from  that  in  a  Jet  flap  system. 

B.  TRANSLATING  FOIL/PARTIAL  CHORD  MECHANICAL  (SUSPENDED  FLAP) 

C.  TRANSLATING  FOIL/PARTIAL  CHORD  FLUID  (SUSPENDED  JET  FLAP) 

Both  schemes  would  use  relative  motion  of  a  part  of  the 
chord  to  reduce  the  excursion  magnitude  of  the  foil.   These 
schemes  seem  worth  further  analysis. 


100 


IV.      LIFT  SPOILERS 

A.  MECHANICAL    (SPOILERS) 

Mechanical  Spoilers  have  been  used  in  aircraft  as 
ailerons.   The  scheme  for  a  hydrofoil  ship  wave  alleviation 
would  require  the  foil  to  be  at  the  angle  of  attack  necessary 
to  produce  the  maximum  lift  needed.  Spoilers  would  then 
reduce  that  lift  down  to  the  minimum  lift  needed.   The  power 
and  weight  for  such  a  system  will  be  small,  but  the  drag  is 
expected  to  be  great. 

B.  FLUID  (VENTILATION) 

The  scheme  of  ventilating  portions  of  the  lift  distribu- 
tion has  been  developed  by  SUPRAMAR.   However,  the  same  pro- 
blems of  drag  enter  as  in  Section  A.   It  may  be  that  overall 
system  power  for  such  a  scheme  would  be  small  compared  to 
full  incidence. 
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APPENDIX  B 

DERIVATION  OF  MOMTBW 
FOR  SECTION  IV-B 
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APPENDIX  C 

DETAILS  OF  THE  IMPLICIT  SOLUTION 
FOR  ALPHA  AND  EEtT  FOR  SECTION  IV-D 
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APPENDIX  D 

PROGRAM  LISTINGS  FOR  SECTION  IV-H 
I.   COMPUTER  PROGRAM  POILDYN  (WRITING) 
II.   COMPUTER  PROGRAM  POILDYN  (PLOTTING) 
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